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Abstract 

In the context of studying black hole singularities by the AdS/CFT 
correspondence, we study the BTZ black hole by a scalar field propagating 
on it and the corresponding two-point Green functions. We explore how 
positions inside the horizon are encoded in the boundary theory. The main 
idea is to relate two different semi-classical approximations of the Green 
function and see how this indicates the bulk-boundary correspondence. 
From a key observation of Festucia and Liu, which is a frequency-geodesic 
identification, we deduce a geodesic approximation from the saddle point 
approximation. As an application, we find saddles of the Green function 
and hence their corresponding geodesies. The conclusion is that some 
of these geodesies do go inside the horizon. This gives the possibility of 
resolving the singularity from the boundary theory. 
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1 Introduction 

Lack of satisfactory description of black hole singularities indicates the limitation 
of classical geometry. String theory as a generalization of general relativity is 
a candidate to resolve the singularities. One way to understand the black hole 
singularities is using the so-called AdS/CFT correspondence conjectured by Mal- 
dacena [T]. This conjecture is motivated by a simple comparison between pairs 
of decoupled theories occurred in some limit of superstring theory. One guiding 
principle of the correspondence is that the identification of isometry group of the 
bulk and the conformal group on the boundary. The proposal is the correspon- 
dence between an isometrically invariant theory on the bulk and a conformally 
invariant theory on the boundary. An extensive review on this subject is given 
in 12]. 

The motivation is simple, while the proof or even the exact statement of the 
conjecture is still far from complete. There are lots of efforts on matching the 
dual theories to support the conjecture. Early examples are [3], HI El [6] . Instead 
of waiting for an ultimate proof, lots of us already enjoy the implication of the 
correspondence. The correspondence in principle allows one to study the bulk 
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theory from the boundary point of view. We will give a brief summary on the 
methods of studying bulk geometry by the boundary CFT. 

The original proposal and its support as [6] are both for the Euclidean sig- 
nature of AdS space. Banks et. al. [5] and Balasubramanian et. al. [7| in- 
dependently give its generalization in the Lorenzian signature. We describe the 
correspondence in the example of a Klein-Gordon free scalar field on the bulk 
AdSd+i- The (positive- frequency) modes solution is decomposed to two classes: 
normalizable modes and non-normalizable modes. On the bulk Klein-Gordon 
quantum field theory, the non-normalizable modes are sources on the background 
and the normalizable modes are those which generate the Hilbert space of physi- 
cal states in the particular background. Based on the principle that dual theories 
should have compatible symmetries required by the identification between isom- 
etry group of the bulk and the conformal group on the boundary, the boundary 
operator related to a KG field on the bulk will be an operator coupled with a 
source on the boundary, which is from the restriction of the non-normalizable 
part of the bulk field. 

At the quantum field theory level, the Hilbert space of the boundary theory 
is generated by the image of normalizable modes from the bulk. The annihila- 
tion and creation operators are the image of those from the QFT on the bulk. 
The corresponding vacua of the dual theories are understood as states that are 
annihilated by all annihilation operators in their respective theories. 

The fundamental relation is the identification between partition functions of 
the dual theories. That is, the partition function of the KG theory on the bulk 
defined by its classical action is equal to the generating functional of boundary 
operator coupled with the source. This identification, in particular, implies the 
correspondence of two-point functions occurring in the dual theories. 

After a general view of the correspondence for AdS space, we now give a 
brief summary on the correspondence for BTZ black hole [8j, which is simply 
a quotient space of AdS^ space. The BTZ/CFT correspondence is studied at 
the level of modes solution of KG-field by Keski-Vakkuri [9] as a generalization 
of [7j. Maldacena in [TU] lays out the BTZ/CFT correspondence in a Hartle- 
Hawking-Israel state. Inspired by Israel's idea in the subject of thermo-field 
dynamics of Schwarzschild black holes [H], Maldacena proposed that an AdS 
Schwarzschild like black hole can holographically be described by considering two 
identical, non-interacting copies of CFT in a particular entangled state. They are 
entangled in the sense that they communicate before their respective evolutions. 
The H-H prescription of gluing half of the Euclidean spacetime to half of the 
Lorenzian spacetime defines the Hartle-Hawking-Israel entangled state. This H- 
H prescription implicitly assumes the analytical continuation procedure, which 
is a defining feature of the original Hartle-Hawking's Feynman propagator on a 
Schwarzschild black hole [12] . Hence it is natural to regard it as a Hartle-Hawking 
like state on a AdS black hole. Under this fixed initial state, one can talk about 
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the correspondence between the dual quantum theories definitely. 

Among various objects in the dual theories, we concentrate on the two-point 
correlation function. In the Euclidean signature AdS/CFT, Klebanov and Wit- 
ten [6] give the sample calculation of obtaining the two-point function (in the 
Poincare vacuum) from the bulk and hence obtain the matching with the two- 
point function on the boundary as a support for the conjecture. The main step 
in Witten's method is to give a bulk-boundary Green's function and so that for 
any boundary field, a unique bulk can be obtained. In passing from Euclidean 
signature to Lorentzian signature, the new feature is appearance of normalizable 
modes [7|. This breaks the uniqueness of finding the bulk field. However, in 
the Hartle-Hawking-Israel state, the gluing process mentioned above implies the 
choice of the bulk-to-boundary Green's function as the analytic continuation from 
that in the Euclidean signature. Therefore, the whole procedure of finding the 
two-point function is just like the procedure in the Euclidean signature without 
ambiguity [THl HU]. The two-point correlation function in the Hartle-Hawking- 
Israel vacuum will be considered as a Hartle-Hawking Green's function of BTZ 
black hole. This is the general framework of the correspondence of two-point 
functions. We will see how to use them to tell the geometry in a holographic 
manner. 

The very first question on regard of studying bulk geometry from boundary 
point of view is how to indicate the radial direction in the boundary theory. Ear- 
lier work in this direction are for example [Ml [15] and more recent attempts are as 
in[l6lll7l[l3l[Tni[lHl[l9l[2Ql|2l]. An important tool in this content is the geodesic 
approximation [17] of the two-point Green function. In the semi-classical limit, 
the two-point correlators can be represented in the form of Feynman paths inte- 
gral of quantum mechanics, and the integration is dominated by certain geodesies. 
In this way, the geometric picture on the bulk appears. It is possible to probe 
singularities by studying the boundary objects which can be related to geodesies 
going inside the horizon. Works on this aspect are for example [T6l [T71 [T9l [20] . 

In the same stream, Festuccia and Liu [22] studied the AdS^ black hole sin- 
gularity. They gave a method to match space-like geodesies in the bulk with 
momentum space Wightman functions of CFT operators in the semi-classical 
limit. 

We study the two-point Green function of a massive free scalar field in the 
Hartle-Hawking-Israel state on the BTZ geometry and its dual on the boundary of 
BTZ. From a bulk point of view, we obtain two-point function essentially by the 
method of modes sum. On one hand, the Green function, written as the Fourier 
integration of the momentum space Green function, can be approximated by the 
methods of steepest descent [23j. The integration is dominated by its evaluation 
at saddles. On the other hand, if we write the Green function in the Feynman 
paths integral of quantum mechanics, it can be approximated by the geodesic 
approximation. The key to relate the two approximations is the observation of a 
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frequency-geodesic relation. After a precise analytic continuation procedure, the 
relation becomes mathematically well-defined. 

This can be explained in the following way. Consider this scalar field prop- 
agating on the BTZ black hole and its two-point function with the two points 
A and B inserted on distinct boundaries in the Penrose diagram. Taking semi- 
classical limit, it can be seen as a wavepacket propagating along geodesies. This 
can also be explained as in the Feynman path integration of quantum mechan- 
ics, where dominant terms are certain geodesies joining A and B. It is plausible 
to make a term by term identification of the WKB approximation of the wave 
equation and the geodesic equation, both appearing as equations of motions in 
the semi-classical limit. After proper analytic continuation procedure, this iden- 
tification gives a definite mapping between frequencies (saddles) and geodesies. 
Hence a geodesic approximation can be deduced from the properties of saddles. 

Mathematically, this frequency-geodesic relation gives the following interest- 
ing property: Firstly, real frequencies parametrize Euclidean time separation of 
geodesies and these geodesies have real turning point outside the horizon. Sec- 
ondly, purely imaginary frequencies parametrize Lorentzian time separation of 
geodesies and these geodesies have real turning point inside the horizon. The po- 
sition space Green function in terms of Fourier transformation has an integration 
contour along real frequencies and saddles picked are near real axis, where the 
corresponding geodesies are with turning points outside the horizon 

Furthermore, if we analogously deform the integration contour from the real 
axis to the imaginary axis we get a new function, assuming that a regularization 
of the function exists. Saddles of the new function correspond to those geodesies 
with real turning point infinitely closed to the singularity. In this sense, we find 
signatures of the singularity on the bulk in the boundary correlation function. 

We should make two remarks regarding the above method: First of all, this 
is by no means the only way to make a geodesic approximation. For example 
different analytic procedure can lead to different dominating geodesies, while the 
correlation function remains the same. This kind of argument was made by Kraus, 
et. al. in |18j. They analytically continued the correlation function directly in 
two different procedures but found the same physical amplitude. Secondly, the 
frequency-geodesic relation does not have physical meaning beyond the semi- 
classical limit. The reason is that the validity of this relation is supported by the 
principle of wave-particle duality. This is only meaningful in the semi-classical 
limit. For example, for a frequency which is not at saddle point, the corresponding 
geodesic may not join the two points in the two-point function. 

The plan for the rest of the paper is as follows. Section two will give a summary 
of the BTZ metric and the some remark on field quantization on BTZ space. In 
section three, we will give the BTZ/CFT correspondence and specify the choice 
of Hartle-Hawking-Israel state and write down a Hartle-Hawking Green function 
by summation of modes. Section four and five will give the explicit mapping 
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between geodesies and saddles on the BTZ black hole at the semi-classical limit by 
comparing the Hartle-Hawking Green function in two different approximations. 
Section six is an application of this geodesic approximation. Finally, we draw 
some conclusions. 

2 The BTZ geometry and field quantization on 
BTZ 

2.1 Classical construction of BTZ black hole 

We summarize the original construction of BTZ black hole of Banados, Hen- 
neaux, Teitelboim and Zanelli [8j and also give it in various coordinate systems. 
Extensive description of the geometry is also in the review [2] . 

The AdS^ space 

Since the BTZ black hole metric is obtained as a quotient metric from the AdS^ 
metric. So we summarize the geometry of AdS^, space first. The AdS^ space is a 
hyperboloid in side the flat R^'^, with metric 

dsl2,2 = -dXl - dXl + dXl + dXl, 

defined by the equation, 

Xl + Xl-Xl-Xl = R\ (1) 

The equation ([T]) can be solved by setting the following coordinate transfor- 
mation from Xj's to (r, p, ^2), 

Xq = R cosh p cos r, 

X3 = R cosh p sin r, 
Xi = i?sinhp cosip, 
X2 = -Rsinh p sin ip, 

where 0<r<27r,p>0 and <■?/'< 27r the solution covers the entire AdSs 
space once. Hence the coordinates (r, p, ip) are called the global coordinates of 
the AdS. The restriction of the metric ds'^2,2 onto this hyperboloid is thus 

ds\as3 = R^ {- cosh^ pdr'^ + dp^ + sinh^c^V'^)- (2) 

This is the AdS^ space. This space has the topology of S*^ x with 5*^ the 
closed time direction in the r direction. To avoid closed time like curve, we may 
unwrap the circle to obtain a straight line R. This space is called the universal 
covering space of AdS, dentoted as CAdSz- 
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Conformal rescaling of CAdS^ as half of ESU 

Since the causal structure has the property of invariant under conformal rescaling 
of a metric, we can study the causal structure of CAdS^ by studying that of a 
conformally rescaled metric, which will turn out to be the Einstein static universe. 

We firstly introduce the following coordinates. We relate to p by tan 6' = 
sinhp, with <6 < n/2. Hence the metric (E]) takes the form 

where — oo < r < cx) and < ip < 27i. This coordinate transformation is in 
particular a conformal transformation of the original CAdS. Now we conformally 
rescaling the metric to a new metric as 

dsl/^ = -dr^ + de'^ + sin^ (4) 

If 9 takes value < < vr rather than < < 7r/2. This is a metric called the 
Einstein static universe (ESU), if the domain of 9 is [0,7r]. Therefore we obtain 
a conformal rescaling of AdS^ as half of an ESU. 

Boundary of the conformal rescaled CAdS^ 

We are interested in the boundary of the conformally rescaled CAdS^ space. In 
the global coordinates of the original CAdS^, the boundary is at p = oo. The 
corresponding boundary of the conformally rescaled (half of the Einstein static 
universe) would be at 6' = 7r/2 in the coordinate system {T,9,'ip). In particular, 
the rescaled metric restricted on the boundary is give by 

ds%/2\e=7r/2 = -dr"^ + dip"^, (5) 

where — oo < r < oo, < ip < 2tt. We see that the boundary is actually a whole 
ESU of dimension two. 

It is not hard to show that the boundary of the conformally rescaled CAdS^ 
as in is actually a conformal compactification of the Minkowski space M^'^. 
It is the identification of the isometry group of the conformally rescaled CAdS^ 
and the conformal group of the conformally compactified Minkowski space that 
plays fundamental role in the AdS/CFT correspondence. 

The BTZ space 

The BTZ black hole metric will be obtained from the quotient metric of CAdS^ 
by certain discrete isometry subgroup of the 5*0(2, 2). In general, if ^ is a Killing 
vector on the spacetime, the transformation from any point P to exp(s^)(P) 
in CAdS^ is an isometry. The group {exp(s^) : s G M} is an isometry group 
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parametrized by s. If s is restricted to s = 27m for n e Z, then it parametrises 
a discrete subgroup. We call it The quotient space CAdS^fEl^^) inherits 

a quotient metric from CAdS^. Under the coordinates (Xq, X^, X^, X3), it is 
proposed that the Killing vector 



up to an isometric transformation. The BTZ metric is defined to be the quotient 
metric induced from CAdSs through the discrete isometry subgroup S(^). How- 
ever, to exclude time-like curves on the quotient metric, ^ • ^ > needs to be 

satisfied everywhere on the spacetime. Hence the final BTZ black hole spacetime 
is defined to be the BTZ metric on the region where C, ■ ^ > 0. These regions can 
be divided into an infinite number of regions of two different types. 

1. Region of type I, < ^ • ^ < 00 (the outer region). 

2. Region of type H, < ^ • ^ < (the inner region). 

We restrict to two continuous regions of both types and introduce the coordinates 
{i, r, -0) as follows. 



1. Region I: r+ < r. 



Xq — — cosh r+ijj, X3 — — — sinh r+i, 
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X — — sinhr+'0, X —— ^coshr+i. 

r+ r+ 



2. Region II: r < r+. 



y> \/ T"^ ~\~ T"^ 

Xq — — coshr+V) ^3 — ~ coshr+t. 



X^ — — sinh r^ijj, X"^ = —— ^sinhr+t. 

r+ r+ 

In the coordinates (f, r, ■0), the Killing vector is ^ = d/dt/j and the discrete sub- 
group acts as taking the identification ^b = ip + 2mr, for n G Z. The BTZ 
black hole metric is thus the quotient metric by the isometry discrete subgroup 
generated by ^ on the region where ^ ■ ^ > 0, 

dslrz = -{r^ - rl)dP + ^^dr^ + r'#' (6) 

r r_,_ 

with the identification ip = ip + 2n7i, n ^ TL. This is a black hole metric, in 
which the black hole singularity is actually a conic singularity, rather than a true 
curvature singularity. 
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Causal structure of BTZ black hole 



To study the global structure of BTZ black hole, we introduce the Kruskal coor- 
dinates as follows, In the region < r < r+. 



U = il ■ smh —t; 

r + r_|_ it"' 



V = \ l ■ cosh —t. 

r + r_|_ K'^ 



In the region r+ < r < oo. 



U = il ■ cosh —if:; 



l^=,/ -^ + ^^ sinh^t. 
In either region, we make the coordinate transformation from {U, V) to {p, q) 

as 

f/ + l^ = tan^^; 

f/-'l/ = tan^^. 

2 

The result is that (i) r = oo is mapped to the lines p = ±|, (ii) the singularity 
r = is mapped to the lines g = ±|, and (iii) the horizon r = r+ is mapped 
to p = ±q. The {U, V) coordinates give the Kruskal diagram and the (p, q) 
coordinates give the Penrose digrams. 



2.2 QFT on BTZ 

QFT in curved spacetime 

To define a quantum field theory in a Minkowski space time, one always has a 
distinct time direction to start with and this time direction together with the 
Poincare symmetry of the spacetime is the key to define a unique quantum field 
theory. In particular, any two inertial observers in Minkowski spacetime agree 
on what is a vacuum state. 

The story on a general curved spacetime is different because it is possible 
to define infinity many inequivalent quantum field theories mathematically |2l] . 
Take the example of quantization a Klein-Gordon field again. If we pick a time 
direction t as in the fiat quantum field theory can be defined as long as 

we find the set of positive frequency modes as eigen modes of the time Killing 
vector and the Klein-Gordon inner product on a Hilbert space spanned by these 
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modes. Vacuum can also be defined by the annihilation operators appeared as 
coefficients of expansion of the field in the basis. However, we should say that 
different choice of t, the respective sets of positive modes may or may not span 
the same Hilbert space. In the later case, we obtain two distinct quantum field 
theories. 

The lack of uniqueness of QFT should not be seen as a drawback curved 
spacetime, rather a sign of the nonzero curvature. The ambiguity arises only due 
to the incomplete (in the sense of unifying gravity and QFT) status of our theory. 
It is also such nonuniqueness leads to the discovery of Hawking radiation of black 
hole [23]. 

Hartle-Hawking state of a Schwarzschild black hole 

Despite the lack of a general quantization procedure of fields on curved space, 
one can still by analyzing particular physical process and make choice of some 
physically meaningful QFT in a curved spacetime. One such example is the QFT 
in Hartle-Hawking state defined on a Schwarzschild black hole geometry. It gives 
the physics of Hawking radiation of black holes. 

There are several equivalent explanation of the Hawking radiation. We would 
like to take a historical routes. Inspired by the nonuniqueness of vacuum states 
in curved spacetime. Hawking constructed a model so that black hole radiating 
particles is possible. For a black hole with its apparent curvature singularity in 
the center, its asymptotic past and asymptotic future are still fiat Minkowski 
space. Therefore, a unique quantum field theory still can be defined on these 
regions. However, the choices of time direction in the region connecting the two 
asymptotic areas are not equivalent, since they cannot be mapped to each other 
by a Poincare transformation. Each choice of time direction will induce some 
ambiguity, which is caused by the curvature. These ambiguities are explained 
as some kind of quantum fluctuations of the curved spacetime. Add them up. 
Hawking proposes the probability of radiation of particles from a black hole |25j . 

Later Hartle and Hawking by generating Feynman paths integral for 
quantum mechanics on fiat space to that on curved spacetime and obtain an 
alternative description of the black hole radiation. The proof is based on the 
construction of a Feynman propagator on the black hole obtained as the analytic 
continuation of its Wick rotated Euclidean Feynman path propagator. The ob- 
servation is that the rate of the black hole radiation obtained in this model is the 
same as the thermal radiation of a black hole located in a thermal bath. 

Inspired by Hartle and Hawking's observation. Gibbons and Perry [26] provide 
another equivalent point of view to Hawking and Hartle's model. Seen from an 
observer standing in distance from the black hole singularity, the problem is 
treated as a grand canonical ensemble of states of the field at a temperature. 
They in particular represent the Hartle-Hawking Feynman propagator in terms 
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of a thermal Green function from the perspective of an outside observer. The 
thermal Green function is obtained by mode sum at Hawking's temperature. 

Israel [Hj provides a third alternative derivation of Hawking radiation based 
on the subject named thermo-field dynamics. In Israel's description, the two- 
point function can be generalized to the case when the two points are inserted on 
distinct boundaries as proposed by Maldacena in [TU] in the course of studying 
AdS black hole in the AdS/CFT. The essential idea is that a quantum field 
theory, with a physical Hilbert space say Ti live on the right side of the Penrose 
diagram can be explained as a composite space 7i ® H where H is defined to 
be a formal Hilbert space live on the left side of the Penrose diagram defined in 
such a way that for any operators from the right side, its expectation value in 
the composite vacuum is the same as its statistical averrage for an ensemble in 
thermal equilibrium. From first sight the dual Hilbert space H is formal, however, 
it can be a full quantum field in its own right. 

Field quantization on AdS space 

Back to the field quantization on curved space time. Since the BTZ black hole 
is a quotient of AdS space, we want to comment on field quantization on AdS 
spaces [271 EHl EHl ED] . There are two obstructions in quantization on AdS space. 
(1) it has closed time-like curve. This can be solved easily by considering CAdS 
space instead; (2) the spacial infinity of (C)AdS is time-like, which means that 
CAdS is not global hyperbolic and hence Cauchy's problem is not well-posed on 
this geometry. It is the time-like boundary that allows information leaking at 
infinity within finite time. This implies that the resulting quantum field theory 
is not a closed system and Cauchy surface can not be well-defined. 

In the case when the field are coupled to the metric conformally, one can trans- 
form the problem of quantization of a field on ESU space, because the CAdS space 
can be conformally rescaled to half of an ESU space. (In the coordinate systems 
we use, the "half" means the range of6'isO<0<7r/2 rather than < < vr.) 
The quantization of field on CAdS can be obtained by the quantization of a field 
live in a box in ESU with the walls at 9 = Tr/2. The boundary condition added on 
the wall could be reflective boundary condition, which will give conservation of 
energy. With this boundary condition added, the Cauchy problem becomes well- 
defined and a unique quantization procedure can be carried out on CAdS. That 
is, the restricted positive-frequency modes which satisfy the particular boundary 
condition span a Hilbert space with the inner product as an integration over the 
Cauchy surface and the quantization procedure is carried out as normal. For 
other possible boundary conditions so to give a well-defined quantization, one 
may refer to references above. 

There are also various quantization procedures provided for AdS black holes 
[311 [321 [331 [311 [351 [361 [371 [38]. Depending on the system, boundary condition on 
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the horizon or origin is needed, too. In our case on the BTZ black hole, we use 
the boundary condition as the reflective boundary condition mentioned above. 
Because it is the normalizable condition in BTZ coordinates which is the 
natural one to impose in the content of BTZ/CFT correspondence. It is also 
the analogue outmodes as used in [2S] in the classical Hartle-Hawking radiation. 
We are not attempting to give a quantization procedure of KG field on BTZ in 
this article. Rather, we specify the particular boundary conditions and obtain a 
Green function at Hawking temperature by summation of modes. This will be 
done explicitly in the next section. The resulting Green function can be regarded 
as Hartle-Hawking Green's function on BTZ black hole. 



3 BTZ / CFT correspondence and Hartle-Hawking 
Green functions 

Assuming that we have a Hartle-Hawking Green function on the BTZ already 
and that the various equivalent explanations of Hartle-Hawking Green's function 
for Schwarzschild black can somehow be carried over to BTZ black hole, we want 
to obtain its various form interested us in the BTZ / CFT correspondence. 



3.1 The BTZ/CFT correspondence 

We summarize the BTZ/CFT correspondence proposed by pUl IT3] . Based on 
Israel's description in the field of thermo-field dynamics on black holes [TT], Mal- 
dacena [10] proposes that the spacetime can be holographically described by two 
identical, non-interacting copies of entangled CFTs on the boundaries. 

Consider the Euclidean metric of the BTZ black hole with metric in the fol- 
lowing form 

dslBTZ = (r^ - ^)dr^ + + r^dcf)^, 

where we have set r_|_ = 1, with the identification = + 2tt, and r = r + /?, 
with (3 the inverse of Hawking temperature of the black hole. Equivalently, if we 
write z = \z\e'^'^, then the metric take the form 

2 . dzdz 

dSEBTZ = 4^Y3^2)2 + 

The boundary of the Euclidean black hole is at r = cxd, the corresponding 
CFT is thus parametrized by (r, 0). It is indeed a CFT at finite temperature 
f3~^. In the {z,T,(f)) coordinate, the boundary is at \z\ = 1. 

The Euclidean signature metric with z = t + r can be analytic continued the 
imaginary part of z as Im{z) = t^, where the imaginary part of z containing the 
Lorentzian time, and now in the Lorenzian section of the whole complex z. If 



(1 + 


\z\ 


2^2 


(1- 


\z 


2^ 
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we write z = —v\ and 2; = m, we obtain the Lorentzian eternal black hole in the 



Kruskal coordinates (u,v,(t)) as 



2 —Adudv (1 — Mf)^ 2 

d^BTZ — T^ \ Y2 ^ 7i \ yi^^ ■ 
(1 + uv)^ (1 + uv)^ 

If we further define coordinates ti.xhy U = t^+x and V = ti — x, so the diagram 
in {tL,x,(f) = constant) is like the usual Kruskal diagram. In this coordinate the 
boundary r = cx) in the AdS coordinates maps to the two hyperbola at UV = — 1, 
the past and future singularity are at UV = 1. 

Back to the complex coordinate z, it contains both the Euclidean section 
and the Lorentzian section. Presumably, the two sections of spacetime are not 
related. Restricting on the Lorentzian signature, we have a particular state at 
time t^ = 0. The two QFTs on the two sections will start to evolve with respect 
to their own Lorentzian time. A particular initial states at tz, = can be simply 
given by the Euclidean section of the space time at Im{z) = (see Fig. 3 of 
|10]). This is referred to as the HH prescription. Thus defined state is called the 
Hurtle- Hawking- Israel state. We will always fix this initial state in our following 
analysis. 

Let us examine what the state is from the boundary point of view. Recall 
that the Euclidean signature AdS^ is a 3-disk. Excluding the S^, it is a 2-disc, 
with the angular direction parametrized by the Euclidean time r. The boundary 
at M = 00 is the boundary of the two-disc (3-disc if considering direction), 
which is the circle. The Euclidean time restricted on constant w-slices always 
have the length jS. The HH-prescription says that the Im{z) = slice in the 
Euclidean black hole is glued to the Lorentzian section 7^/2 x at each constant 
M-slices. When restricted this identification at the boundary, (in the Euclidean 
signature the boundary is m = cxd half of the circle surrounding the two-disk; 
in the Lorentzian signature the boundary is simply the hyperbola UV = — 1), 
we obtain two Lorenzian boundaries connected by half of a Euclidean circle. 
The proposal is that two CFTs live on two pieces of the Lorenzian boundaries, 
with the initial states equivalent to a state live on half of the Euclidean circle 
I 13/2 X S^. These are two decoupled, entangled CFTs. They are entangled from 
their communication at the start. 

If we denote the two copies of the Hilbert space of the two CFT by Hi and 0.2, 
and the whole Hilbert space Ti = Tii x I-L2 of the system. The initial wavefunction 
1^) G 7i of the Lorentzian evolution is of the explicit form 

= ^L=Y,^M-PEj2)\Er,)i X |E„)2, (7) 

where the sum runs over all energy eigenstates and Z{(3) is the partition function 
of one copy of the CFTs at temperature I3~^. ([7]) is the Harlte-Hawking-Israel 
state. 
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3.2 Hartle-Hawking Green functions 



Hartle-Hawking Green function in the bulk and the boundary 

With the specification of the particular vacuum states, we consider, in the semi- 
classical approximation, the correlation functions for insertions of operators on 
distinct boundaries [TOl |T3]. For the BTZ black hole the correlation functions 
of points on distinct boundaries are decided by the analytic continuation, whose 
validity is guaranteed by the particular choice of the Hartle-Hawking- Israel state. 

Suppose a two-point Hartle-Hawking Green function with two points on the 
same boundary, say the one on the right, is given as (^((r, (r, 
(r — > oo) where = -|- t and m_ = — t. The analytic continuation implies 
that when one point moves to the left boundary, the new Green function is 
obtained by substitute entities in G by 

i(3 

u± — >u±T y • 

It is this analytic continuation implicitly commits one to consider a particular 
states of CFT defined on the disconnected BTZ boundaries and it is implied by 
the HH-prescription mentioned before. 

Similar procedure can be carried out from the boundary point of view, we 
will denote the Green function in the boundary theory by G. The simple relation 
between the bulk and boundary Hartle-Hawking Green functions is 

G{{u+,u.),{u'u'^)) = lim z-^''+z'-^''+G{{z,u+,u.),{z',u'u'_)), (8) 

z,z' >0 

where = 1=*=^!+™-^ is the conformal dimension of the field restricted on the 
boundary. 

Following the discussion in Section 12.21 we will give an explicit two-point 
Green function for the bulk by summation of modes. Due to the boundary 
condition of modes we choose, this can be understood as a Hartle-Hawking two- 
point function on the bulk. 

Modes solution of the Klein-Gordon equation 

In the following, we will find the normalizable modes mentioned above for a 
scalar field on the BTZ black hole. The Klein-Gordon equation of a free scalar 
field <l>(r, t,0) propagating on the BTZ metric ([U]) is given as 

□$ = m^<l>, (9) 

where □ is the Laplacian operator defined by the metric. We can solve this 
equation by the method of separation of variables. Write the ansatz of positive 
frequency modes as 

fuj,jit,r,(j)) = exp(iJ0) exp(-iu;it)r~^/^X^,j(r), 
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where J is integer. 

The radial part of the equation of motion is 



- - 1)'^%^ - 2Kr^ - 1) - uj'X^Ar) + Vir)X^Ar) = 0, (10) 
where 

V{r) = (r^ - 1) + "'-1)^ v' = l+m\ 

Related to terms in Section [3l2l we have /i± = (1 ± z^)/2. 

A standard quantization procedure follows this would be (1) add sensible 
boundary condition on these modes at infinity and only pick a subspace of these 
positive frequency modes to construct a Hilbert space; (2) specify the Klein- 
Gordon inner product on the Hilbert space and further normalize the basis to 
obtain an orthonormal basis; (3) expand the field solution in this basis and pro- 
mote the corresponding coefficients to annihilation and creation operators; (4) 
determining vacuum state by using annihilation operators; (5) the corresponding 
Green function in this vacuum can be obtained by summation of these orthonor- 
mal modes. Unfortunately, it will be beyond the scope of this paper to finish 
these steps rigorously. To obtain the final Hartle-Hawking Green function we 
will only involve (1) and (5). Such construction is by no means unique. 

Equation ffTUjl can be found in terms of hypergeometric functions, 

X^,j(r) =r-2{l- r^y'f {Cir-'-^G{r) + C2r'-^H{r)) , (11) 

where 

. „ / -iuj - iJ - u + 1 -iuj - iJ + u + 1 ^ . ^ ^ 
G{r) = F , ; 1 - zJ, 



H{r) = F i , ; 1 + tJ, 

and Ci,C2 are constants. 

We specify the boundary condition for the positive- frequency modes as follows. 



Xiz) 



exi){iujz) exp(i(5) + exp{—iujz) exp{—i5) as z — > oo (horizon) 
C{uj, J)z^^'^ as z — > (boundary) 

(12) 

where S is some real constant and the coordinate z is defined by 

z{r) = - In (r + l)/(r — 1). 
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The boundary condition at infinity agrees with the reflective boundary condition 
in [39] in global coordinates. It is also the same as the boundary condition 
satisfied by normalizable modes in [9]. Therefore, from the point of view of a 
well-defined QFT and an AdS/CFT related modes, this boundary condition is 
natural. Furthermore, these modes when restricted to U = {V = 0) have an 
expansion in terms of e~*'^^(e~*'^^) with u > and hence are related to Kruskal 
modes in the same way as out modes related to Kruskal modes, which is the 
defining modes of the Hartle-Hawking state of a Schwarzschild black hole. With 
all these inspections, we may regard the summation of these modes at Hawking 
temperature as a Hartle-Hawking like Green function of the BTZ black hole. The 
boundary condition we imposed at the horizon seems less transparent, except 
technically easy, and needs a justification of its physical meaning. 

With these boundary conditions, we are able to fix constants Ci and C2. The 
detailed calculation is in the appendix. The coefficient 

„, -,,0 exp(7rcc;)(exp(27rco') — 1)^^^ /I , . 

r ( -(-iuj + iJ + u + l)]r (-(iuj -iJ + 



2' 7 V2 

T (^^{-iuj - iJ + u + 1)^ (13) 

at boundarjQ, which will be used. 

Hartle-Hawking Green function by modes sum 

The Wightman function at zero temperature is given by the modes sum. 



G+((t,r,0),(t',r',0')) 



/oo 
c/cu^exp(^J(0' - 0) - tiu{t' - t))r-^/V-^/2 
■00 J 

Note that u is not quantized under the boundary condition (fT2|) . On the other 
hand, the Fourier integration representation of G^{{t,r,(f)), (t',r',0')) is 

/oo 
dtu J2 exp(a(0' - 0) - iuj{t' - t))g+{uj, J; r, r'). 
■00 J 

Comparing coefficients in the two integrations above, the Wightman function in 
the momentum space is 

^+(^, J; r, r') = r-'/'r'-'/'X^j{r)X^j{r'). (14) 



^ This constant is worked out with the help of Mukund Rangamani. 
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The Wightman function Q'^{uj, J;r,r') at the Hawking temperature Th = ^, 
where /3 = 27r, can be obtained from the Wightman function at the zero temper- 
ature by [ID], 



1 — exp{—(3uj) exp{l3uj) — 1 



This is the momentum space Hartle-Hawking Green function. We can pass the 
relation between position space two-point functions (Ej) to the momentum space 
two-point functions. The boundary Green function can be obtained as, 



^+(u;,J) 



r,r' >oo 

hm r^+'^r' 

' >oo 

exp(/3co') 



exp(/3u;) 



r,r' >oo exp{pu) — 1 



exp{/3uj) 



-c{iu,jy 



Substituting (1131) in, we obtained the boundary thermal Green function in the 
momentum space, 

;;4-/ exp(37rti;)ti; ^ /I , ^ n\ ^ /I / 



(15) 



ly 



u=x+iy 



Figure 1: Poles of Qt{u, k) in the large u limit. 
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To analytically continue this function by extending the domain of uj from the 
real line to the complex tu-plane, we need to identify singular points of ^^(t^, J\ 
They are located at 

{cu = ±J ± (1 - v)i T 2M} U {cu = ± J T (1 - v)i ± 2M}, (16) 

where N is any positive integer. For later use, we change the scale of variables as 
u = uj/v and k = J/v and further assume that A; is a positive real number. This 
assumption will be made clear later by relating it to the assumption of L = iLi 
where L/ > in Section 14.11 Singular points fll6p in the large v limit are located 
at 

{u = ±A; ± z ± 2m} U {m = ±A; =Fi =F 2m}, (17) 

where n = N/ v and is any positive integer. 

On the complex w-plane, the distance between poles become very small as 
V — > oo and the four lines of poles become branch cuts. This is illustrated in 
Fig. [H These are the branch cuts for analytical continuation of Q^[u^ k). 

4 Geodesic approximation and Saddle point ap- 
proximation of the Green function I 

We want to see how the exact boundary correlation function probe the bulk 
geometry. The main idea is to relate correlation functions to geodesies on the 
bulk. To achieve this, we consider two different semi-classical approximations of 
the two-point Green function. The first method is the geodesic approximation. 
Indeed, wavepackets propagate like particles along geodesies in the semi-classical 
limit. 

With respect to the decomposition of complex time separation into the Lorentzian 
and Euclidean sections, i.e. t + it, the BTZ spacetime can be projected onto the 
two sections as shown in Fig. [21 

We are interested in the case when the two points are inserted on distinct 
boundaries of the Penrose diagram, describing the Lorentzian section of the space- 
time. Geodesies joining distinct boundaries are like A, B, C and D in Fig. [2] 
(a). 

The second approach is by the steepest descent method to evaluate the cor- 
relation function as Fourier integration of the momentum space Green function. 
For points inserted on distinct boundaries in the Lorentzian section space-time, 
the boundary two-point function can be obtained from the Wightman functions 
by changing the time variable from t' to t' — i(3/2. That is, 

Gu{t,t'-k)= [ exp(-tuuU' -t-tf^jjgf3+{u,k). (18) 



oo 
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(a) 



(b) 



Figure 2: (a) Space-like geodesies joining points on distinct boundaries with 
time separation t' — t — ij3 /2 in the Lorentzian section of the spacetime. Lines 
of equal Lorentzian time are straight lines through the centre, the time being 
given by the slope. Lines of equal radius r are hyperbolas in each region (not 
shown), with asymptotics at the cross in the middle. The upper and lower waved 
lines are the future and past singularities. The upper part of the cross is the 
future horizon r = 1 and the lower part of the cross is the past horizon. The 
left and right straight lines are two distinct boundaries at r = cxd. Each wedge 
is associated with a constant Euclidean time. Region I is of time 0, region II is 
of time — i/9/4, region III is of time — and region IV is of time —32/9/4. The 
jumping of — ?/9/4 comes from the fact that integration of proper time from one 
region to the next picks up a pole at the horizon with a factor of —if3/4. A, B, C 
and D all have Euclidean time separation —(5/2. As u varies from —ioo to 
along the imaginary axis, geodesic moves from Aio B and then to C. Geodesic 
D corresponds to m/ > 0. (b) Space-like geodesies in the Euclidean section. The 
radial direction of the spacetime is along the radial direction of the circle. The 
centre of the circle is chosen to be at the horizon, r = 1 and the boundary of 
circle is where r = oo. The Euclidean time r is in the angular direction with 
period 2ti. As u varies from +oo to 0^, the corresponding geodesic changes from 
A to -B, and then to C . Note that C turns back before it reaches the horizon at 
the centre. Geodesic D stands for a geodesic with m < 0. 
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The asymptotic expansion of the integration is dominated by saddle points. The 
goal is to relate the saddles and the geodesies, both appeared as dominating 
objects in the semi-classical limit, so that we can deduce a geodesic approximation 
from the saddle point approximation. In other words, we want to make the 
following geodesic approximation explicitly defined in the semi-classical limit, 

Gu{t,t'; k) = 5^exp(-m(l, + /,(A;))), (19) 

i 

where i indicates different dominant geodesies, m is the mass of the scalar field 
and Ij is proper length of the i-th dominant geodesic joining the two points and 
fi{k) is a function indicating the contribution from the angular momentum k. We 
will firstly give the notion of geodesies on the BTZ black hole, secondly sketch 
the method of steepest descent and thirdly provide a formal proof of f|T9|) . We 
will deduce an explicit geodesic approximation in the next section. 

4.1 Some formula related to geodesies on the BTZ black 
hole 

On the BTZ metric ([6]), there are two particular Killing vectors, dt = (1,0,0) 
and dtf, = (0,0, 1). Their corresponding conserved quantities are energy E"^ and 
angular momentum L^. Since energy and angular momentum are conserved on 
geodesies, we can write down geodesic equations in terms of E and L. Let 
x(A) = (t(A), r(A), 0(A)) be a space-like geodesic with the affine parameter A. Its 
geodesic equation can be written as 

1 fdry E^ 

where E = {r"^ — l)dt/d\ and L = r^d(l)/dX. If we define an effective potential by 

Veff = -(r^ - 1) (l - ^) , (21) 

then ( l20i) gives {dr/dX^ = E"^ — Vg//. This can be thought of as the equation of 
motion of a particle of energy E"^ in a potential K//- Observe that the potential 
(l^Tj) behaves differently according to the different values of L. 

See Fig. [3l Imagine a particle of energy E'^ coming from boundary towards 
the singularity. When L is a real number, the potential is represented by the 
Curve 1 of Fig. [3l We can see that with big enough E"^, the particle will go 
inside the horizon and be trapped there. However, when L is a purely imaginary 
number, the potential is represented by the Curve 2. We can see that no matter 
how big E"^ is, the particle has to bounce back. When E'^ is infinitely large, the 
turning point of the geodesic becomes infinitely closed to the singularity. In this 
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1 1 ii>+T|rHH-|44444ij_rri 

t 1.0 1.5 2.0 2.5 3.0~ 



Curve 1 
Curve 2 



Figure 3: Potential v.s. radius with different L. 



case, tlie returned particle may carry signatures of the singularity. Therefore, we 
want to consider the case when L is a purely imaginary number. For simplicity, 
we assume L = iLj for Lj > 0. The turning point of such particles can be 
solved as solutions of the equation 



(22) 



The proper time t, proper length I and proper angular displacement d of a space- 
like geodesic with energy E"^ and angular momentum are as follows. 



t{E,L) = 2 



E 



(r^-l)JE^ + ir^-l){l-^) 



\{E,L) = 2 f 

POO 

d(E,L) = 2 / 



1 



+ (r2 - 1) (1 - ^) 
L 

E2 + (r2 - 1) (1 - ^) 



■.dr. 



(23) 



(24) 



(25) 



where the integrations are all over the contour along the real r axis from to 
positive infinity. 
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4.2 Steepest descent methods 

We will briefly summarize the steepest descent methods from [23l HI]. It is a 
technique to get asymptotic expansions at the large v for integrations in the 
form, 

^W) — / giu) exp{iyh{u))du, (26) 
Jc 

where C is a contour on the complex w-plane, and g{u) and h{u) are holomorphic 
functions. The points of w-plane where dh{u)/du = are called saddle points. 
They are saddle points on the real surface representing | exp{i'h{u))\, which is seen 
as a function of the two real variables x and y with x = Re{u), y = Im{u). Curves 
along which Im{v{h{u)) is constant are called steepest descent paths. Along such 
curves, | exp(z//i(-u))| changes as rapidly as possible. In other words, they are 
the gradient lines of | exp(z//i(M))|. The saddle point is a stationary point of the 
function \ ex\i{yh{u))\ along the steepest descent path. The method of steepest 
descents consists in deforming the path of integration C so as to make it coincide 
as far as possible with arcs of steepest paths. In this way, the integration is 
transformed to a real integration of the Laplacian type and the Laplace method 
may be used to evaluate the integral asymptotically. We summarize the methods 
in the following steps: 

1. Identify the saddle points Uq as zeros of dh{u)/du = 0. 

2. Determine degree n of saddle point uq such that 

dih(u), d^'hiu), , , 

l"o = 0^ g = 1,2, . . . ,ra-l, = aexp(za), a>0 

for some real number a. 

3. Determine directions of steepest descent of each saddle point: \i u — Uq = 
pexp{i9), then steepest descent directions are 

a (2n+l)7r 

e = — + ^ ^, p = 0, 1, . . . , n - 1, 

n n 

and steepest asscent directions are 

e = + -^, p = 0,l,...,n-l. 

n n 

4. Justify, via Cauchy's integral theorem, the deformation of the original con- 
tour of integration C onto one or more of the paths of steepest descent. 

5. Determine the asymptotic expansions of the deformed integrals through 
Laplace methods. 

For our interests of finding dominating saddles and eventually their corresponding 
geodesies, we do not need to go through the fifth step. The rest of the steps listed 
will become clearer when we consider our particular examples. 
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4.3 Saddle points and the geodesic approximation 

We will formally deduce ( |T9l) for now and provide a mathematically explicit 
description in the next subsection. Consider the Fourier integration of the two- 
point function (1181) in the large v limit, 

Gi2{t,t'; k) = J duexp ^—iuu ^At — i^^^ exp{2h'Z{u, k)), (27) 

where At = t'—t and we write Q^{u, k) = exp(2z/Z(ti, k)) for some Z{u, k) so that 
the Fourier integration is transformed into the form of We can analytically 
continue the function Q^{u, k) from the real w-axis to the complex w-plane with- 
out difficulty by branch cutting pole lines illustrated in Fig. [H Complex saddles 
can thus be obtained by solving the saddle point equation. The evaluation of the 
integrand at the dominating saddles is dominant in the asymptotic expansion of 
[T8l On the other hand, we ca use the language of Feynman paths integration 
of quantum mechanics for (1271) in the semi-classical limit. There will be certain 
geodesies joining the two points to dominate the integration in the spirit of wave- 
particle duality. Since the saddles and geodesies appear respectively as dominant 
terms in two different semi-classical approximations of the same correlation func- 
tions, we want to make a link between them. The first attempt is to simply make 
a term by term identification, by comparing units, between the wave equation 
and the geodesic equation [22] ■ 

To obtain the wave equation in the semi-classical limit, we use the ansatz 
X{z) = exp{i>S{z)) for the radial part of the Klein-Gordon field. Equation ( fTOl) 
can thus be written in terms of S{z). We substitute in the expansion S{z) = 
So{z) + l/uSii^z) + l/u'^S2{z) -|- ■ ■ • , in which the Sq{z) will be the leading term 
in the large v limit. The radial part of the Klein-Gordon equation can further be 
written in terms of 5*0 (z). 



dSoiz) 



--, + -y^u' = 1. (2^ 



r"^ — 1 \ dz J x 

On the other hand, recall the geodesic equation (I2UI) . describing space-like geodesies 
joining points on the boundary of energy E"^ and angular momentum L^. The 
assumption of equivalence between equations (l28l) and (!20|) suggests the following 
identification: 

dSo dr 

— — < — > — , k < — > iL, u < — > lE. (29) 
dz dX 

Observe that the relation between L and k is consistent with the assumptions 
that L = iLj for Lj > and that k is positive. 

Now we will show that fl^Ul) implies the formalism of the geodesic approxi- 
mation of propagators as in ( |T9l) . Write the momentum space Green function in 
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terms of So{z), i.e., 



Ggiu,k) 



lim 



exp{j3uj) 
exp(/3a;) — 1 



r ^/^r' exp(z/ 5*0(2;) ) exp (z/ 5*0(2:') )■ 



(30) 

and define ZQ{u,k) = lim^(r) ,qSq{z). The semi-classical approximation of mo- 
mentum space Green function is thus Q'^{u, k) ~ exp(2z/Zo(u, k)) where u » 1. 
We apply the identification fpUl) and change coordinates from 2; to r in the inte- 
gration. We obtain 



Zo{u,k) = lim / 



lim 




r2 - 1) 1 + — - m2 



111 
= iu- -t{u,k) - -\{u,k) + -L(i{u,k), 



where r^. is the solution of 



+ 



^2 ^ 



(31) 



(32) 



and t, I, d are from fl23|) . fl24|) and fl25|) . respectively. Also notice that equation 
agrees with equation (1^ under the assumption (j^HD- Here "=" reminds us 
that these are formal equalities. We will fix an analytic continuation procedure 
to make these equalities exact later. We make some immediate remark regarding 
(I3T1) . It implies that 



2dZo{u,k) 
d{—iu) 



-t{u,k), 



2dZo{u,k) 
d{-ik) 



d{u, k). 



If we regard the proper length l(t, d) as a function of (t, d) instead of {u, k), then 
the above identities imply that 2ZQ{u,k) and l(t, d) are related by a Legendre 
transformation. Geometrically, a Legendre transformation maps the graph of a 
function to the family of tangents to the graph. In our case, if the graph is l(t, d) 
then the tangents are c?l(t, d)/dt = iu. Alternatively, if the graph is 2Zo(u, /c) 
then the tangents are 2(iZo(M, k)/du = —it. As Lagrangian and Hamiltonian we 
have the mathematical relation 



2Zo{u, k) + l(t, d) = iut - ikd 



and hence 



d2Zo{u,k) 
du 



it. 



(33) 
(34) 
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We further consider the decomposition of the identification —iu = (r^ — l)dt/dX 
into the real and imaginary parts, and take integration respectively. We obtain 

tE = J ^P^d\ + CE, tL = j -^d\ + CL, 

where Ce and are constants of integration. We will fix them by fixing a 
reference geodesic Gq. We choose it as the geodesic with time separation — z/5/2, 
which is the same as geodesic C in Fig. [2] (a). So knowing ur and ui is equivalent 
to knowing Ie and II- In other words, we assign m = to the geodesic Go and 
use other u to measure how other geodesies are different from Gq in proper time 
separation. For example, if = 0, then the change of the proper time from that 
of Go is zero in the Lorentzian section. If ur is equal to zero, then the change of 
the proper time from that of Gq is zero in the Euclidean section, that is to say 
the Euclidean time separation will always be the constant —i(3/2. 

Now come back to the semi-classical approximation of Q^{u,k) in (1301) . By 
applying flHTl) . we get 

Q~p{u,k) ~ exp(z/(m ■ t('u. A;) — !(«, A;) + Ld('U, A;))). 

The position space boundary two-point Green function ( |T8l) is therefore, 

Gi2(t,t';k) ~ J du exp ^—iu — t — 

exp (z/ {iu ■ t{u, k) — \{u, k) + Ld{u, k))) . (35) 



Evaluating by steepest descent methods, the integration fl35|) will be dominated 
by its evaluation at saddle points, which are solutions of 2dZo{u, k)/du = i{At — 
iP/2). By identity ( IMl) . the saddle point equation implies 

t = At-i^. (36) 



That is to say, saddle points as frequencies under the identification (1291) will be 
identified as geodesies with proper time separation equals to At — i[3/2. This 
agrees with the physical picture we have in mind. However, under the same 
identification, frequencies which are not saddles correspond to geodesies with 
time separation not necessarily At — i/3/2. This means that those geodesies may 
not join the two points which appeared in the two-point function. The reason 
is that the relation ( 129|) is obtained in the semi-classical limit when the wave- 
particle duality holds. We have no reason to expect it to be physically meaningful 
beyond this limit. Finally, by substituting into and applying (IHUj) . we 
have 

Gi2{t,t';k) ~ J (i-u exp ^— i-u ^At — z/^ exp(z/(— l(t, d) -|- iwt — i/cd)) 
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where the last term is obtained as a summation with respect to saddle points, 
labelled by z, and identification between v and m is assumed. The previous 
fi{k) which appeared in (fT9|) is obtained as ikdi. Hence we induce a geodesic 
approximation (fT9l) from the assumption of the identification (129|) . 

5 Geodesic approximation and Saddle point ap- 
proximation of the Green function II 

Since ( l29l) is a complex relation, the correspondence we obtained in the last sub- 
section is only a formal correspondence unless a proper analytic continuation 
process is provided. We have already seen that the Green function can be an- 
alytically extended by specifying its branch cuts. We will provide an analytic 
continuation for the geodesic approximation and hence make the formula fll9l) 
rigorous. There are two steps involved to define an analytic continuation for the 
functions t, I, d. The first step is to analytically continue rc{u) to an analytic 
function on the complex w-plane, and the second step is to deform the contour 
of integration to the complex r-plane. 

For the first step, is a solution of flH2|) . or equivalently, r{uY + (/c^ — v? — 
l)r{uY — k^ = 0. If the discriminant (/c^ — u"^ — 1)^ + 4k'^ = 0, then u is in the set 

{u = ±k±i}U {u = ±kT i}- 

Observe that these points are located at the four end points of branch cuts of 
momentum space Green function at the large u limit as in Fig. [H In this case, 
r{u) has a root of multiplicity 4. While the discriminant is not zero, there are 
four distinct roots. We will only concentrate on the one given as, 

re(u) = (-A;2 + u2 + l + /(A;2-n2-l)2^4p). (37) 

It can be analytically continued to the complex w-plane if we cut the points 
where the discriminant is zero. Since properties we are concerned with do not 
differ much qualitatively when k varies as a positive real number, we will fix A; = 1 
from now on. We plot the real part and the imaginary part of the analytically 
continued function of (137|l in Fig. HI 

We want to mention some properties of the analytic function rc{u). When u 
approaches imaginary infinity, the norm of the turning point |/'c('w)| approaches 
and hence so does rc{u). As u approaches real infinity, the norm of the turning 
point approaches positive infinity. As u approaches from any direction, the 
norm |rc| approaches one. u = corresponds exactly to our reference geodesic 
Go- III particular, as u approaches in the direction along either the real axis 
or the imaginary axis, the imaginary part of rc{u) vanishes. Hence the analytic 
function rc{u) approaches the horizon at r = 1 in these two cases. 
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(a) (b) 

Figure 4: Fig. (a) is the real part of r^lu) on the complex u = x + iy-plane when 
k = 1. Fig. (b) is the imaginary part of rc{u) on the complex w-plane when 
k = l. 



To understand the function rc{u) mathematically is quite simple. As the way 
we obtained, it is a holomorphic function coming from an analytic continuation of 
some function on the real line. Certain limiting behaviours mentioned above are 
quite neat. Meanwhile its physical explanation is not as clear as its mathematical 
explanation. The first problem which occurs is that the concept of turning point 
loses its usual physical meaning. For example, with a particular turning point of 
a complex value, we don't know how to compare it with the horizon at r = 1. 
Take its norm? Take its real part or anything else? So when we see this analytic 
function rc{u) we keep in mind that it is a nice function coming from a real 
function, which has physical meaning. If we want to use it, we use it in a way 
related to the original real function. Now consider the complex w-plane, the 
complex r-plane and the analytic function Tc which maps u to r. Physically, 
there are two distinct positions on the r-plane: the point r = 1 is the event 
horizon and the origin r = is the singularity of the black hole. The previous 
observation concerning the function rc{u) tells us that it approaches r = 1 when 
u approaches zero along any direction. Therefore, we consider simply all straight 
lines, u = \u\ exp{i6) with some constant 6, through the origin on the w-plane 
and find their images under rc{u) on the complex r-plane. 

Particular cases of this mapping are plotted in Fig. O 

Consider points on the line u = \u\exp{i6) starting from zero to the com- 
plex infinity. The analytic function rc{u) thus gives us a family of geodesies 
parametrized by their turning point rc{u) as u moves. No matter what 6 is, the 
starting geodesic is the one with turning point at the horizon, which is Go- There- 
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(a) 



(b) 



Figure 5: (a) shows straight hnes u = \u\exp{i6) with varies slopes, (b) shows 
the images of rc{u) of the hnes in (a). Lines and their images are of the same 
colour. 

fore, we find a way to put any geodesic with an arbitrarily complex rc{u) into a 
family of geodesies by joining the origin to the u. This family is considered to 
consist of deformation of the reference geodesic Go- When ^ = 0, m is always real, 
and the family of geodesies are physically observable geodesies with real turning 
points starting from the horizon (Gq) to positive spatial infinity. When 6 = 7r/2, 
the line parametrizes a family of geodesies with all real valued turning points, 
starting from horizon {Go)and decrease monotonically to the spatial singularity 
at r = 0. Incidentally any line with 6 greater than 7r/4 (this particular angle 
actually depend on our assumption of = 1, other k will give other angle) has 
turning point rc(u) = as m goes to infinity. (Fig. [5l) 

As to the second step, the original contour in the integrations t, I, d as in fl2^ . 
([21]), fl^ starts from rc{u) and goes to infinity along the real r-axis. As rdu) is 
analytically continued into the complex r-plane, we need to specify the contour 
of integration. Although the deformation of the contour can be arbitrary, we 
would like to make a choice so that the heuristic relation ( |T9l) gives a proper 
geodesic approximation of the correlator. Since when rc{u) is very far from the 
horizon, the choice of contour is not easy to make, we will start by considering 
when rc{u) gets close to the horizon at r = 1. Write rc{u) = l + eexp(z?7) where e 
is a sufficient small positive number and rj ranges from to 27r. Recall that rc{u) 
and M = |m| exp(i^) are in one-one correspondence as shown in Fig. Hlor Fig. [5l 
We will decide our contours according to the way how rc{u) approach to 1 with 
respect to rj and hence to 6. As u approaches to 0^ from +oo along the real 
axis {6 = 0), they are proposed to relate to geodesies with difference in proper 
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time from the geodesic Go only in the Euchdean section. Their turning points 
vary from r = +00 to r = 1. For example, geodesies A, B and C plotted in (b) 
of Fig. O This is actually the case before analytical continuation, so we simply 
choose the original contour as in (a) of Fig. [61 In this case, we can see that the 
pole at horizon does not contribute. This immediately implies the choice of the 
contour for u approaches to 0~ from —00 along the real axis {9 = vr) as in Fig. [H] 
(b). Indeed, simply consider a closed contour in (a) of Fig. [21 it starts at r = 0, 
r = 00, and then goes very close to the horizon and finally goes to r = 00 ( as C). 
After this, it moves downwards along the arc a little bit and comes back along 
the dashed geodesic, which is very close to the horizon. Finally it goes back to 
the point where it starts. Integration along this closed contour will give proper 
time difference by 27ri. While we have already chosen the contour as (a) of Fig. 
[6l so that the outwards geodesic avoids the singularity, we have to make a choice 
so that the backwards geodesic can pick up the factor of 2ni. The contour in (b) 
of Fig. [HI when circling around the pole, picks up this factor. Hence the choice 
is confirmed. 

That leaves us two more directions to consider. For the purpose of making a 
geodesic approximation, we are interested in geodesies joining two boundaries as 
A, B, C{Go) and D in the Penrose diagram in Fig. [21 These geodesies correspond 
to u = iuj, since there are no Euclidean time difference between them and the 
reference geodesic Gq. The imaginary part uj is indeed related to the Lorentzian 
proper time of the geodesies, since the Lorentzian time separation of Go is zero. 
We want to make a choice of contour so that the jumping in the Euclidean time 
— i/3/2 is included even if we only carry out the integration for the Lorentzian 
section. As Uj moves from +ioo to along the imaginary axis, rc{ui) moves from 
r = to the horizon along the real axis. As rc{ui) goes very closed to r = 1, we 
choose the contour to be (c) of Fig. [HI so that this naturally include the jumping 
of — i/9/2 by making use of the pole. Symmetrically, as uj moves from —ioo to 
iO~ along the imaginary axis, the contour is chosen as (d) of Fig. [6l 

Therefore, for the turning point very closed to the horizon, we have made our 
choices for 9 = 0, 7r/2, tt, 37r/2 as in (a), (c), (b), (d) in Fig. [HI respectively. We 
make the following requirement for the rest of the 6''s: when < 9 < 7r/2, the 
contour is chosen to be as in (a); when 7r/2 < 9 < n, the contour is chosen to 
be as in (c); when it < 9 < 37r/2, the contour is chosen to be as in (b); when 
3it/2 < 9 < 2tt, the contour is chosen to be as in (d). As to u away from 0, i.e. 
\u\ not very small, rc{u) is far from the horizon. We simply choose the contour 
according to the argument 9 oi u. 

Combining the two steps, we obtain a complete analytic continuation pro- 
cedure. In this way the previous heuristic relation (JT9l) is exactly defined. A 
interesting phenomenon is that the Lorentzian section of geodesies, which is re- 
lated to the imaginary part of the saddle points, can have turning point inside 
horizon. In this sense, we may be able to find a signature of the BTZ singu- 
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Figure 6: Choices of deformed contours, (a) The contour chosen for u approaches 
to O"*" along the real axis, (b) The contour chosen for u approaches to 0~ along 
the real axis, (c) The contour chosen for u approaches to iO"^ along the imaginary 
axis, (d) The contour chosen for u approaches to iO~ along the imaginary axis. 
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larity from the boundary correlation function. However, we by no means imply 
this is the only possible geodesic approximation. It is possible to assign different 
correspondence for u and Vc firstly, then to define a corresponding analytic con- 
tinuation according to its relation with the correlation functions secondly. That 
is to say, it is possible to define a different way of geodesic approximation for 
the same correlation function. Similar ideas can be seen in [TH], in which two 
different but equivalent analytic continuations for position space Green function 
are presented. 



6 Application of the frequency-geodesic relation 

With the exact geodesic approximation provided in the last section, we would like 
to see how it works through a simple computation and its physical implication. 
We will first look for saddle points of the correlation function fl7r|) and then look 
for their corresponding geodesies. There are two cases of interest. The first case 
is when the time separation is purely Euclidean as in (b) of Fig. O The second 
case is when the two points sit on different boundaries in the Lorentzian section 
of space time with time separation At — if3/2 where At is the Lorentzian time 
separation as in (a) of Fig. [2l 

In the large u limit, instead of finding So{z) and hence Zq{u, k) as we did be- 
fore, we adopt a slightly different procedure. Since we already have the explicit 
expression of the momentum space correlation function, we obtain the approx- 
imation simply by considering its asymptotic form in the large u limit. After 
rescaling u, J to u, k, we can write (|T8|) in the form of fl27|) and get the expression 
of 2Z{u, k) = l/h'\nQi3'^{u, k) as 

+ -\n\T \^\iu-ik + l+ 



Applying the Stirling formula |42j to the Gamma functions appearing in 2Z(m, k\ 
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we get 

= I (■" { 2.v(l+ + ^'"^^ - - ^ + (i)) + 

+3™ + -(i-u + i/c + 1) ln(i-u + iA; + 1) 

+ — ik + 1) \\i{iu — ik + 1) + -^[—iu + ik + 1) ln(— m + ik + 1) 
+-(— m — ^/c + 1) ln(— m — ik + 1). 

The derivative with respect to u and is then obtained as follows, 

2dZ(u,k) i, (iu + ik + l)(iu — ik + 1) , , 

yLJ. = 37r + -ln-^ — '—, asz/»1.38 

du 2 {-tu + tk + l){-tu - tk + 1) ^ ' 

To carry out the steepest descent method, we may write (J2711 in the form of ( l26l) . 
In this way, giu) will be simply the constant v and 

= -iu (^t - i0 + 2Z(m, A;). 

Saddles of 6*12 (At — ijS k) are thus determined by dh{u)/du = 0, that is, 

2rfZ , /3 

^ = iAt + - 39 
du 2 ^ ' 

We write At = t + iT for t and r both real numbers, i.e. t is the time separation 
in Lorentzian section and r is the time separation in the Euclidean section. We 
can see that t = is our first case where the time separation is purely imaginary, 
r = corresponds to the second case. From ( l38l) and ( 139|) . we can find saddle 
points by solving the following equation, 

{iu + ik + l){iu — ik + 1) 



[—iu + ik + 1)(— m — ik + 1) 



exp(2(t + ir)). 



where the inverse Hawking temperature j3 = 2tt is substituted in. There are two 
solutions, 

, i{l + exp(2At)) ± v/-4exp(2At) + k^{l - exp(2At))2 
= 1 - exp(2At) • 

For simplicity, we will again assume k = 1. Computation can be carried out for 
other non-zero positive k similarly. When t = 0, we can view saddles as being 
parametrized by the Euclidean time r. Since the Euclidean time separation r is 
of period 2tt, we may assume r G (— vr, tt]. The saddles are 



g ^(l + exp(2zr))± v/-4exp(2zr) + (l-exp(2ir))2 

■u, r = ; ^ . 40 

^ l-exp(2ir) ^ ^ 

32 



We observe that (t) have vanishing imaginary part. This imphes that the 
time separation is purely Euchdean, the possible dominating geodesies of the 
correlation function will not have Lorentzian time section. Another observation 
from the solutions is that the two saddles are both periodic functions with respect 
to r of period vr. In particular, u^{t — tt) = u^{t). So there is in fact only one 
solution as seen in Fig. [71 




Figure 7: parametrised by Euclidean time r. As r varies from —tt to 0, the 
graph of u^ij) is shown. As r varies from to tt, the graph of u^ij) is shown. 

Besides the formula ( HOl) . the behaviour of the saddle can be seen clearly from 
the graph of the correlation itself as in the left column of Fig. [HI As the Euclidean 
time difference r increases from to tt, the saddles move from to the right along 
the real w-axis. For the limiting case when r = tt, the saddles actually move to 
positive infinity. The changing of the shape of the graph is with respect to the 
moving of the saddle. Symmetrically, when r varies from to — tt, the saddle 
starts from the origin and move to left till negative infinity along the real axis 
(not shown in the figure). 

By evaluating the second derivative dPh{u)/dv? at the saddles, we find that 
the saddle in this case is of degree one for all r G (— vr, vr] and the steepest descent 
contour is the real axis itself. Thus in this case, there is no need to deform the 
integration contour, and the saddle obtained dominates the integration. 

According to the frequency-geodesic relation specified in the last section, we 
can find the corresponding geodesies with the saddle points. These geodesies will 
play the role of dominating geodesies in the geodesic approximation of the two- 
point function. As illustrated in Fig. [9l as m increases from 0, the corresponding 
geodesic moves from H'H, to B' B and then A' A. As the saddle increase, the 
turning point of the corresponding geodesic move from the horizon to infinity. 
This is also what we expect. 

The second case we want to consider is r = 0. Saddles are parametrized by 
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(e) (f) 

Figure 8: Graphs of |exp(/i(u))| on the complex u-plane. (a), (c) and (e) are 
parametrized by the Euchdean time separation r when t = 0. (a) is when r = 0, 
the saddle is at m = (0,0). (c) is when r = 1.5. (e) is when r = 2.5. (b), (d) 
and (f) are parametrized by the Lorentzian time separation t when r = 0. (b) 
is when t = 0. (d) is when t = 0.1, the upper saddle moves downwards and the 
lower one moves upwards rapidly, (f) is when t = 1.5, the two saddles stops at 
the lines x = ±1. the lower part of the graph reduce height rapidly while the 
upper half increases its height rapidly. 
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H' 



Figure 9: As saddle moves from the origin to the positive real axis, the corre- 
sponding geodesic moves from H'H{i.e. Go) to B'B and then A' A. 



the Lorentzian time t and given as 



L ^ ^(1 + exp(2t)) ± v/-4 exp(2t) + (1 - exp(2t))2 

' l-exp(2t) ■ ^ ^ 

Saddles are plotted as in the Figure [TOl For each time t between and oo 
there are two saddle points. As t varies from to around 0.88, the upper saddle 
moves from the origin downwards to —i and the lower saddle moves from negative 
imaginary infinity upwards to —2i rapidly. The two saddles "meet" when the time 
separation is about 1.10. After that, the upper saddle moves to the left and the 
lower one moves to the right, both with constant imaginary part. After the time 
separation exceeds t = 1.7, the two saddles disappear into the branch cuts around 
u = ±1 — i. 

A magnified picture indicating the "collision" of the upper and lower saddles 
is shown in Fig. [TTl The behaviour of saddles can also be seen from the graph of 
the correlation function itself as shown in the right column of Fig. [HI When the 
Lorentzian time separation is zero, the upper saddle is at the origin and the lower 
one is at —ioo. There are two "lumps" (both diverge at infinity) in the graph. 
We refer them as the upper and lower part. As the time separation increase, the 
lower part descend (this corresponds to the lower saddle, who always sits on top 
the "lump" , moving upwards) and the upper parts at the same time grow. When 
the time separation becomes sufficiently large, the lower part almost becomes fiat 
and the upper part grows very tall. 

These are both of the saddles solved as solutions of dh(u)/du = 0. However, 
not all saddle points will contribute in the final approximation. We need to find 
out their steepest descent direction and decide which one can be used in the 
approximation. If we evaluate the second derivative d^h{u) / dv? at the saddles 
obtained, we can determine the steepest descent direction for each saddle. This 
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Figure 10: Two trajectories of the upper and lower saddles u^{t) parametrized 
by Lorentzian time t (increases in the vertical direction) on the complex w-plane. 

is also indicated in Fig. [HI We can see that before the two saddles collide, the 
upper saddle has steepest descent direction in the horizontal direction. The lower 
saddle has its steepest descent in the vertical direction. Since the contour of the 
Fourier integration is along the real w-axis, we conclude that the dominant saddle 
is the upper one. In the cases when the saddles are colliding and after collision, 
their positions and steepest contours are symmetric with respect to the y axis, 
so both of the saddles contribute. 

Now we would like to see how movements of saddle points imply the geodesic 
approximation. At t = 0, the upper saddle is at the origin u = 0, and this 
corresponds to the reference geodesic Gq. This geodesic is shown as HH' in (a) 
of Fig. [121 The lower saddle is at imaginary infinity and does not count as a 
saddle yet. When the Lorentzian time separation t increases (from H'H, to A' A, 
B'B and finally to C'C in (a) of Fig. [T^ . the upper saddle moves downwards; 
the lower saddle moves upwards. Correspondingly, for solid geodesies (related to 
the upper saddle, which is dominant), their turning points move from the horizon 
to inside the horizon. For dashed geodesies (related to the lower saddle, which 
is not dominant), their turning points come out from the singularity but remain 
inside the horizon. When the Lorentzian time separation increases to C'C, which 
is related to the "collision" of the two saddles, the solid geodesic and the dashed 
one almost coincide. Consider projections of these geodesies onto the Euclidean 
section. Since the relating saddles both have vanishing real part, both projections 
will have time separation exactly — 1/3/2 as the reference geodesic Gq. Therefore, 
for any geodesic appearing in (a), we can plot the centre of of the circle to be 
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(a) 



(b) 




(e) (f) 

Figure 11: Saddles move with respect to Lorentzian time separation t. (Suggest 
to read the left column first and then the right one.) Shading in the figure has 
the property: darker means lower and lighter means higher, (a) t = 0.60. (c) 
t = 0.80. (e) t = 0.88. (b) t = 0.90. (d) t = 1.10. (f) t = 1.70. The upper saddle 
dominates in (a), (c) and (e). Both saddles dominate in (d). Saddles disappear 
on to the branch cut in (f). 
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Figure 12: Geodesies in the Lorentzian and Euclidean sections before the saddle 
points collide, (a) The Lorentzian section. Geodesies with solid lines correspond 
to the upper dominating saddle at various time separation. Geodesies with dashed 
lines correspond to the lower not- dominating saddle at various time separation, 
(b) The Euclidean section at time separation t. When t = 0, X'X = H'H and 
Tc = 1. As t increases, X'X stands for A' A, B'B and finally CC, with their 
corresponding turning points. 

its turning point Vc- Then the Euclidean projection of the geodesic will be X'X. 
X'X stands for H'H, A' A, B'B or C'C, and Vc depends on geodesies specified. 

Before the collision of the saddles, the dominating geodesic is always the solid 
one. After the saddles collide, the saddles pick up real part frequencies at the 
same time and both saddles make an equal contribution in the final approxima- 
tion. (See Fig. [lU]). Correspondingly, there appear two dominating geodesies. 
They both have Euclidean time difference with respect to the reference geodesic 
Go- However, since their differences is of opposite sign, the summation of their 
Euclidean time difference with respect to Gq remains 0. For both of the geodesies, 
their turning points are no longer real valued. So we need to extend both the 
Euclidean section and the Lorentzian section in a third dimension to indicate the 
direction of the imaginary part of the turning points as shown in (a), (b) of 
Fig. [131 The magnified version of the third dimension is shown in (c) of Fig. 
[T^ When the saddles collide, it is indicated as C'C in (a) of Fig. O The two 
very close geodesies still lie inside the two-dimension Penrose diagram, with their 
real turning points inside the horizon. Let solid geodesies correspond to the left 
saddle, which is from the previously upper one and dashed geodesies correspond 
to the right saddle, which is from the previously lower one. As time separation 
increases from C'C, to E'E, F'F and then to C'C, turning points of the solid 
geodesic now slowly rise above the Penrose diagram. Their turning points have 
almost constant real part and increasing imaginary part. This is shown as the 
upper half of (c) in Fig. [131 At the same time, the turning points of the dashed 
geodesic now slowly descends to underneath the Penrose diagram. They are again 
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Figure 13: The spacetime related to the cases of colhsion and after colhsion of 
the two saddles, (a) The Lorentzian section. Solid eodesics correspond to the 
left saddle, which is from the previous upper one. Dashed geodesies correspond 
to the right saddle, which is from the previous lower one. (b) The Euclidean 
section. After the collision of saddles, there are two geodesies dominating si- 
multaneously. The sohd geodesic corresponds to the left saddle and the dashed 
geodesic corresponds to the right saddle, (c) After collision, the turning points 
pick up imaginary part of opposite sign. Real part of is in the horizontal di- 
rection and imaginary part of is in the vertical direction. As time increases, 
the two points start from the real axis, and split to upper and lower direction in 
a symmetrical way. This is the precise image of the vertical arrowed lines in (a) 
and (b). 
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with almost constant real part, but decreasing imaginary part. It is shown as the 
lower half of (c) in Fig. [131 

The Euclidean sections of these geodesies are less simple to describe. For the 
pair of dominating geodesies, (except for C'C,) in (a) of Fig. [131 there are two 
corresponding geodesies in the Euclidean section as their projections. The two 
images have opposite time differences with respect to —i^ and with same real 
parts and opposite imaginary parts of turning points. They are indicated in (b) 
of Fig. [131 The solid geodesic is the projection of the solid geodesies from (a) 
and the dashed geodesic is the projection of dashed one from (a). The centre of 
the disk is decided by the real part of the turning points of various geodesies. 

In summary, for the time separation being Euclidean as ir — i/3/2, one saddle 
point appears. This implies that for the geodesic approximation, there is one 
dominant geodesic for each r. This geodesic has purely Euclidean proper time 
ir — ij3/2. As r increase, the turning point of the geodesies varies from the 
horizon to the infinity. For the time separation being t — if3/2, there are always 
two saddle points appearing as solutions of the saddle point equation. However, 
they don't both dominate all the time. As t increases from to the time when 
the two saddles collide, only the upper saddle dominates. The corresponding 
geodesic approximation tells us that the dominating geodesic starts from the one 
with time separation — 1(3/2 and turning point at the event horizon, i.e. Gq. As 
t increases, the turning point of the geodesic moves inside to the horizon. At 
the same time, the turning point of the other (not dominant) geodesic is moving 
from deep inside the singularity outwards. We will make more comment in the 
next paragraph on the non-dominating partner. By the time the saddles collide, 
the two geodesies almost coincide. After that, two geodesies contribute equally 
in the geodesic approximation. They always have the same real part of turning 
point and opposite sign imaginary parts of the turning point. As t increases, 
real parts of their turning points stay inside the horizon not changing much, 
around 0.7 to 0.8 (see (c) of Fig. [T3l) . while imaginary parts increase in same 
amount in norm. In the Euclidean section, the two dominating geodesies have 
time difference —i(3/2 + i6 and —i(3/2 — i6, rather than before when there is only 
one dominating geodesic with time separation —i(3/2. We have seen that some of 
the dominating geodesies depending on the time separation have turning points 
inside the horizon. In this sense, we say that the boundary correlation function 
encodes locations inside the event horizon. On the other hand, in all the cases, 
the turning points never go infinitely closed to the singularity. This surely is not 
good enough for the purpose of resolving singularity. 

As we mentioned, when the time difference is t — ij3/2 for t quite small, there 
appears the lower saddle. It moves from —zoo upwards along the imaginary 
axis. Its steepest descent direction is always along the imaginary axis. Recall 
that at the same time the upper saddle always has steepest descent direction in 
the horizontal direction. This implies that when the steepest descent method is 
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applied, the lower one is not picked up because of the contour of integration is 
along the real axis. We can define a function from the original Fourier integration 
by replacing the contour of integration by the imaginary axis as an analogue of the 
new observable proposed in p2] for AdS^ black hole. Despite the fact that the new 
integration actually blows up and assuming a regulization is available, we study it 
gularized mathematical function. If we apply the steepest descent method 
as before, the lower saddle rather than the upper one will become dominant. That 
is to say, if this new function has any physical meaning to anybody, wavepackets 
can move along the geodesic with turning point deep inside the horizon at the 
semi-classical limit. The smaller the time separation t is, the deeper the particle 
can go. As t approaches zero, the turning point approaches infinitely close to the 
singularity. In this way, the newly defined function as a function on the boundary 
carries information from the singularity and may help in resolving the BTZ black 
hole singularity. At the same time, we should admit that its physical meaning is 
not clear. 

7 Conclusion 

The main objective of this paper is to see how the two-point function in the 
boundary CFT carries signatures from inside the horizon of the BTZ black hole 
in the context of the AdS/CFT correspondence. By using the key identification 
between frequencies and geodesies in the semi-classical limit in [22], we deduce a 
geodesic approximation from the saddle point approximation. As an application, 
we compute the Hartle-Hawking Green function by mode summation directly and 
obtain the boundary two-point function by the AdS/CFT correspondence. By 
taking the semi-classical limit, we find saddle points of the two-point function 
with respect to the time separation and thus relate the saddle points to specific 
dominating geodesies on the bulk. There are two interesting features which ap- 
pear. The first occurs when the time separation is t — i/?/2, for the Lorentzian 
time separation t not very small. There are two dominant geodesies with turning 
points complex conjugate of each other and also with Euclidean time separations 
— i/?/2 -|- iS and — i/9/2 — i6, respectively. This indicates the approximation of 
the two-point function is given by two complex geodesies together. The second 
interesting feature appears when the time separation is t — if3/2 for the Lorentzian 
time separation t not very far from zero. If we replace the original integration 
contour in the Fourier integration of the position space two-point function by 
the imaginary axis and regularize the new integration, then the saddle picked up 
with respect to the new contour is the one who can have turning point infinitely 
closed to the singularity. If there is good explanation about the newly defined 
function in the boundary theory, it may help to resolve the singularity. 
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Appendix 

W will fix constants Ci and C2 so that the boundary conditions (fT2|) are satisfied. 
We will get an approximate form of the solution around boundary first and fix 
the constraint between Ci and C2 by imposing the boundary condition (|T^ . As 
r — > 00, z ^ — > 0, the approximated solution is 

Xbdryiz) = {-l)-fz-lz'^{C^z''G{z)+C2Z-''H{z)), 

where 

G{.z) = F{ , ; 1 - tJ, — ), 

iJ-iuj + l-u iJ-iuj + l + v 1 
H{.z) = F{ , ; 1 + iJ, —). 

From general identities of hypergeometric function ^42j , we can further approxi- 
mate the solution by 

XMry{z) = (-l)-^((Ci(-l)-^ri + C2(-l)-'^r2)^^^ (42) 

+(Ci(-i)-^r3 + C2(-i)-^r4)^^+'^), (43) 



where 



r4 



r(i -ij)r(z/) 

r(i(i -tuj-ij + z/))r(i(i + iuj-tj + z/)) ' 

r(i + a)r(z/) 

r(i(i -tiu + ij + z/))r(i(i + iuj + tj + u)) ' 

T{l-iJ)T{-iy) 

T{1{1 - lio - tJ - u))T{l{l + iuj - tJ - u))' 

T{l + iJ)T{-u) 

r(i(i -iuj + tj- z/))r(i(i + iuj + tj- u)) ■ 
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Imposing the boundary condition by requiring Z2 ^ mode of ( H2l) vanish, we 
get the constraint between the constant C\ and C2 as 



C2 = -i-^y^C,. (44) 

Thus, ( l42i) becomes 

Xhdry{z)=C{uJ,J)z"^+\ (45) 

where 

C{u;,J) := Ci(-i)^(r3-^r4) 

J- 2 

= Ci(— 1)^^— sin(7rz/) sinh(7rJ) 

Ti-u)T + + r(l - ^J)T ^ ^ + ^- + ^J + - ^ . 

Next, we want to get the approximation of ffTTj) around the horizon and impose 
boundary conditions to fix constants Ci and C2 completely. 

Again by general properties of hypergeometric functions, the approximation 
of ffTTj) around the horizon is 



where 



r(i -ij)r(zcj) 


r(i(i + 


iu — iJ + z/))r(|(l + 

r(i + ij)r(icu) 




^)) 


r(i(i + 


iuj + zJ + z/))r(i(l + zcj 

r(i - zj)r(-zcu) 




^)) 


r(i(i- 


iuj ~ iJ — z/))r(|(l — icij 
r(l + iJ)r(-ia;) 


- iJ + 


^)) 


r(|(i- 


ioj + iJ — z/))r(|(l — zcj 




^)) 



r8 = 

In the z coordinates, 

Xhor{z) = exp(-^)[(Cir5 + C2r6)2-^-exp(-«a;^) 
+ (Cir7 + CaFg) exp(7ro;)2^'^ exp(iwz)]. 
By the constraint ( l44l) from boundary, 

X,„,(^) = exp(-^)[(Cir5-(-l)*V^ir6)2-^'^exp(-za;^) 

+ (Cir7 - (-l)^-'^Cir8)exp(7r^)2^-exp(2a;^)] 
J- 2 

= Ciexp(-^)(2-''^Pexp(-icu2) +2^'^Qexp(icjz)), 
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where 

(oxp(Ja ) - oxp(-Ja ))r(l - tJ)T(t^') 



(exp( Jtt) + exp(7r(a; + ii/)))r(|(l + iu; - iJ - i/))r(|(l + iu;-iJ + u)) ' 

exp(— 7rc<j)(exp( Jvr) — exp(— J7r))r(l — iJ)r{—iuj) 



^ (exp(j7r) + exp(7r(— + iz/)))r(i(l — ito — iJ — z/))r(i(l — iuj — iJ + u)) 

Numerically, it turns out that \P\ = \Q\. The boundary condition at the horizon 
implies that Ci should be chosen as Ci = exp (^)(Pg)-2. Indeed, 

X„or{z) = ^^exp(-^)(2--Pexp(-ia;z) + 2-Qexp(^a;z)) 

[q 

' — exp{—iu}z) + 2*'^y — exp{iu}z). 

If we define 6 by exp(i6') = "J^, then \P\ = \Q\ implies that 6 is real. Let S = 

a;ln2 + ^, then the solution is Xhor{z) = exp{iujz) exp{iS) + exp{—iujz) exp{—iS), 
and hence satisfies the boundary condition. After Ci is fixed, C{u!,J) is fixed 
and its square is given as, 

C{u,JY = -^exp(7ra;)(exp(27ra;) - l)a;sin^(7ri/)r(-i/)^ 

1 1 

T{-{iuj + iJ + U + l))T{-{-iuj + iJ + U + 

Zi Zi 

V<X{iu -iJ + u + l))r{-{-iu} -iJ + u + 1)). 
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